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ABSTRACT 

Many of the most luminous extragalactic radio sources are located at the centers of X-ray clusters, and 
so their radiation must be scattered by the surrounding hot gas. We show that radio observations of 
the highly-polarized scattered radiation (which depends on the electron density distribution) in combi- 
nation with the thermal Sunyaev-Zcldovich effect (which measures the electron pressure distribution), 
can be used to determine the radial profile of the electron temperature within the host cluster. The 
sensitivity levels expected from current instruments will allow radio measurements of mass-weighted 
cluster temperature profiles to better than a ~ 1 keV accuracy, as long as the central radio source 
is steady over several million years. Variable or beamed sources will leave observable signatures in 
the scattered emission. For clusters with a central point source brighter than ~ 1 mJy, the scattered 
polarization signal is stronger than competing effects due to the cosmic microwave background. 
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1. INTRODUCTION 

Massive galaxy clusters often host a bright radio source 
at their center (e.g., Zhao, Burns & Owen 1989). A small 
fraction of the radiation emitted by the central source is 
expected to scatter off the surrounding intracluster elec- 
trons and obtain a high level of polarization. The radial 
profile of the scattered radiation depends only on the 
electron density distribution within the cluster, in dif- 
ference from the thermal Sunyaev-Zeldovich (SZ) effect 
which depends on the electron pressure distribution (see 
Carlstrom, Holder & Reese 2002 for a recent review). In 
this Letter, we show that the temperature profile of the 
intracluster gas can therefore be determined by means of 
radio observations alone, through the combined detection 
of the polarized scattered light from a central source and 
the SZ effect. This method for inferring cluster temper- 
atures is particularly powerful for distant clusters, where 
the alternative option of using an X-ray telescope re- 
quires very long integration times and is less practical. 

The typical optical depth for Thomson scattering 
through the core of mas sive X-ray clusters is ~ 0.01 
(|Mason and Mvers ll2000|) . and so roughly a percent of 
the observed flux from the central radio source should 
appear in the form of a diffuse radio halo around it. To 
give a concrete example, a point source observed to have 
a flux density of I Jy in a cluster at a redshift ~ 0.2 
would generate a ~ 10 mJy radio halo on an angular 
scale of ~ 1', assuming it is emitting isotropically. De- 
spite its faintness, the diffuse radio halo could be sepa- 
rated from contaminating sources due to the high ampli- 
tude (~ 50%) and the lar ge-scale tangential alignment 
of its polarizat i on pattern J Sazonov and jj unyaev 1999; 
iSunvaev lll98l iSunvaev and Zeldovich Ill980|) . The typ- 
ical intrinsic polarization of individual point sources is 
only a few percent, making the central point source a 
much less significant contaminant in the polarization 
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map than it is in the intensity map. With current in- 
struments such as t he VLA achie ving stability over a dy- 
namic range of 10 5 l)Perlevlll999|) . the halo signal should 
be detectable. 

Several groups are currently constructing radio tele- 
scopes optimized to search for the SZ signature of clusters 
of galaxies ( e.g., Carlstrom, Holder & Reese 2003) . The 
planned SZ surveys will be efficient at finding clusters 
of galaxies up to relatively high redshifts [z ~ 1), from 
which the received X-ray flux would be weak. In addition 
to the reduction in the X-ray flux due to the large cosmo- 
logical distance, the prominence of galactic outflows (due 
to star formation and quasar activity) at z ~ 1 — 2 will 
likely disperse somewhat the i ntracluster gas and make 
the cluster les s X-ray luminous l|Evrard and Henry 111993: 
iKaiser 1119911) . Thes e problems are not nearly as acute 
for SZ observations ijHolder and Carlstrom | |200D. rais- 
ing the possibility of strong SZ sources that are difficult 
to follow up with X-ray telescopes. Our proposed method 
provides an important alternative route to probing the 
density and temperature distributions of the intracluster 
gas in the most massive clusters that will be discovered 
by future SZ surveys. 

Most SZ-detected clusters will probably not have a 
suitable central point source for obtaining detailed im- 
ages of the reflected images. In nearby clusters, roughly 
one quarter of clusters host a radio source with radio 
power greater than 10 24 - 5 W/Hz at a wavelength of 21 cm 
l)Ledlow and Owen1ll995() . At z ~ 0.5 this corresponds 
to roughly 10 mJy, so it might be expected that the po- 
larized emission might be detectable in a sizable frac- 
tion. However, the obser ved radio luminosity function 
l)Ledlow and Owe"n"lll99fil) at high power falls off quickly 
(dN/dlnPoc P~ 15 ), so the fraction of sources with high 
enough radio powers to give a strong signal is likely a few 
percent, depending on how much the cluster radio source 
luminosity function evolves with redshift. From recent 
observations at 30 GHz (see Bennett et al. 2003 for a 
recent compilation) of radio point sources, one would ex- 
pect roughly 300 sources brighter than 0.1 Jy per stera- 
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dian. Dee p SZ surveys expect rou ghly 30000 sources per 
steradian ijCarlstrom et al. Il2002j) . so if every bright ra- 
dio source is embedded in a cluster the reflected emission 
should be significant for roughly 1% of detected SZ clus- 
ters. 

The projected profile of the scattered radio emission 
is a direct probe of the electron density run with radius 
n e (r). By comparison, the X-ray emissivity is sensitive 
mainly to n 2 , while the SZ effect probes n e T, where T(r) 
is the radial profile of the electron temperature. The 
combination of SZ measurements with radio polariza- 
tion measurements allows a direct determination of the 
electron temperature profile, while the combination of 
X-ray observations with polarization halo measurements 
allows a direct distance determination, based on the same 
method that is currently us ed with X-ray and SZ obser- 
vations (R eese et aT"ll2002j) . Finally, the combination of 
X-ray, SZ, and polarization observations would allow a 
consistency check on the underlying simplifying assump- 
tions, such as the spherical geometry of the cluster or the 
steadiness of the central source. 

Our proposed method could be complicated by tem- 
poral variability of the radio source, because of the ge- 
ometric time delay that the scattered light acquires. 
Images of extended radio jets indicate that extended 
radio sources are often stable over t hose timescales 
(jBeeelman. Blandford. and Rees I Il984j) . However, if 
central cluster sources do vary on short timescales, then 
an independent mapping of the electron distribution 
around these sources (e.g. based on deep X-ray observa- 
tions) could be used to read off their temporal variability 
as a function of age from the brightness distribution as a 
function of projected radius in their scattered radio sur- 
face brightness. If the source is strongly beamed in a 
direction not aligned with our line of sight, it is possible 
to use the reflected emission as a probe of both the source 
luminosity and the surrounding medium (see, for exam- 
ple Chiaberge et al. 2003 for an application using X-ray 
and ultraviolet imaging). However, for simplicity, we fo- 
cus in this publication on the subject of steady sources 
emitting isotropically. 

In §2 we calculate the polarization of the scattered 
radio halo, while § 3 outlines a synthetic observational 
procedure to test the prospects for temperature measure- 
ments by the above method. We close with a discussion 
of some observational and theoretical hurdles which must 
be overcome as well as a summary of the many poten- 
tially important applications of measuring diffuse radio 
polarization in galaxy clusters. 

2. SCATTERED RADIO HALO AROUND A CENTRAL POINT 
SOURCE 

We start with the idealized case of a steady point 
source of observed spectral flux density S v and spectral 
luminosity L v (per observed frequency v) at the center of 
a spherical cluster. We define a coordinate system cen- 
tered on the cluster with the z-axis along the line of sight 
and calculate the circularly symmetric brightness profile 
of the scattered radiation. An infinitesimal cell of physi- 
cal area dxdy on the sky contains dN e — n e dxdydz elec- 
trons and therefore has a physical effective cross-section 
to incoming photons of dN e aT, where ot is the Thom- 
son cross section. The total spectral luminosity that will 
therefore be scattered by this cell is dN e (aT/^Trr 2 )L I/ , 



where r is the 3D radius from the cluster center. The 
radiation scattered along the z-axis toward the observer 
will be preferentially polarized in the direction perpen- 
dicular (on the sky) to the cylindrical radius vector. This 
is particularly easy to see for a parcel of electrons lo- 
cated in the plane of the sky but offset from the central 
source. An initially unpolarized source will have equal 
components of the electric field E oriented in the plane 
of the sky (perpendicular to the radius vector) and along 
the z-axis. Scattering of the z-polarized photons into 
our line of sight will be suppressed, leading to the scat- 
tered radiation being completely polarized in the direc- 
tion perpendicular to the radius vector in the plane of 
the sky. In general, the suppression factor of the compo- 
nent along the cylindrical radius vector is proportional to 
cos 2 9, where 6 is the angle between the spherical radius 
vector and the z-axis. 

Defining an azimuthal angle on the sky <f> and a co- 
ordinate system for the Stokes parameter Q (in units of 
intensity, not flux) such that the maximum of Q is ob- 
tained at (f> = 0, we find (see also Sazonov & Sunyaev 
1999, Sunyaev 1983) that for a general distribution of 
electron number density n e {r) 

Q{R^) = ^(?f)cos24> f + °°dz^fn e (R^,z) 3 -l 
16?r \JyJ J-oo r sr 

where R is the cylindrical physical (not angular) radius 
in the plane of the sky, r = V R 2 + z 2 , and dA(z) is 
the angular diameter distance to the source (arising here 
from the conversion of physical area on the sky, dxdy, to 
solid angle). 

For an electron number density profile following 

the "/3-profile" , where n„ = n eQ (l + r 2 /^) 3 / 3 / 2 

ijCavaliere and Fusco-Femiano I Il976|) . one obtains for 
the Stokes Q parameter [see, Gradshteyn & Ryzhik 1980, 
3.259(3), p. 299] 

Qifr *) = 15? ( § ) cos 2 <t>W (2) 

x 2 F 1 [2,i;|/? + 2;l-^^] 

xB[I,|(l + /3)] * 

where B(x,y) is the Beta function and 2-F1 (a, b;c;z) is 
the hypergeometric function. A similar expression in 
terms of a Beta function and 2-F1 can be obtained for 
the total scattered intensity. For the case of j3 = 3/2 
an explicit expression can be found in Sunyaev (1983). 
Note that the source redshift does not enter explicitly 
into this expression, since the intrinsic source luminosity 
and reflected luminosity share the same redshift factors. 

Typical parameter values are dA ~ l/i _1 Gpc for a 
source at a cosmological distance, and n e ~ 0.01 cm~ 3 , 
r c ~ 150/i -1 kpc for typical cluster cores. For a 1 
Jy source, these parameters correspond to roughly 1 
mjy/arcmin 2 at about a core radius (~ 0.5') from the 
central source. At an observing frequency of v = 30 
GHz, the relevant conversion to brightness temperature 
is 1 Jy/arcmin 2 = 0.4 K, and so the scattered intensity 
is equivalent to nearly a mK polarized temperature sig- 
nal. In a signal-to-noise sense, this is an easy signal to 
detect for the current generation of instruments. The 
difficulty is that such a signal arises from a field that by 
construction has a 1 Jy source at the center which must 
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be stably removed. Typically, the 1 Jy source may also 
intrinsically produce tens of mJy of polarized flux. Note 
that the polarized scattered radiation will be more cen- 
trally concentrated than the projected gas distribution 
due to the 1 /r 2 fall off of the source flux as seen by the 
cluster electrons. 

3. PROSPECTS FOR TEMPERATURE MEASUREMENTS 

We examine the feasibility of temperature measure- 
ments by constructing synthetic observations of galaxy 
clusters with steady point sources embedded at the cen- 
ter and emitting isotropically. We first use /3-models and 
then examine clusters in three-dimensional cosmological 
simulations with hydrodynamics. Throughout our dis- 
cussion, we consider a 1 Jy source observed at 30 GHz. 
All properties of the polarized radio halo, scale directly 
with the point source flux. The scattered intensity from 
the point source may exceed the SZ signal amplitude in 
the core of the cluster. In principle, high-frequency SZ 
observations arc not necessary, since the scattered in- 
tensity can be reconstructed from the polarization and 
removed from the SZ flux. 

First we assume a /3-profile for the radial profiles of the 
gas temperature and density, and construct the bright- 
ness profiles for the SZ effect and polarized radio halo 
around the central source. Rather than tailor our predic- 
tions closely to a particular instrument, we assume that 
instruments such as SZA 2 , CARMA 3 , AMI 4 , Amiba 5 , 
EVLA 6 , or ALMA 7 may generically be thought of hav- 
ing sensitivity to a range of angular scales from roughly 
arcsecond scales up to a few arcminutes with up to 
sensitivity. Note that SZA and AMI are not currently 
planned to be used for polarization studies, but it would 
not be difficult to adapt these instruments for polariza- 
tion studies. 

The details will depend crucially on the observing fre- 
quency. In general, large scale emission will be filtered 
out. We assume that this filtering is on a scale of ten 
arcminutes or a radius in the u-v plane corresponding 
to ~ 300 wavenumbers (in units of the observing wave- 
length) and that there is no sensitivity to larger scales 
(i.e., we assume a high pass filter). At cm to mm wave- 
lengths this is effectively required to avoid contamina- 
tion from primary cosmic microwave background (CMB) 
anisotropies. In addition, high angular resolution is re- 
quired to isolate the central point source. We assume a 
maximum u-v radius of 3000 wavenumbers and an equal 
noise contribution per log interval in u — v radius be- 
tween the above inner and outer cutoffs. Noise was added 
such that the resulting noise in each map was 1 fiK and 
beam effects were ignored. An observing frequency in 
the Ray leigh- Jeans was assumed for the SZ component. 

We assumed a /3-profile for each of the electron num- 
ber density profile and the pressure profile, with a central 
electron number density of 0.01 cm~ 3 , a central temper- 
ature of T e0 = 5 keV, a common core radius of 0.75' for 
both n e {r) and T(r), [3 e = 2/3 for the electron density 
profile and /3p = 1 for the pressure profile. The cluster 

2 http:/ /astro. uchicago.edu/SZE/survcy.html 

3 http://www.mmarray.org 

4 http:/ /www.mrao.cam.ac.uk/telescopes/ami 

5 http:/ /www.asiaa.sinica.cdu.tw/amiba 

6 http://www.aoc.nrao.edu/evla 

7 http://www.alma.nrao.edu/ 
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Fig. 1. — Reconstructed temperature profile for a single realiza- 
tion (assuming the electron density and temperature profiles are 
both /3-models) with approximate l-cr uncertainties. Neighboring 
points are highly correlated. The true input profile is shown as a 
solid line. 



was assumed to be at z — 0.5. A realization of "data" 
was generated and then fit to a model with two /3-models, 
with the centers of the profiles constrained to be at the 
cluster center. Each /3-model had three parameters, and 
the presence of a point source at the center was allowed 
in terms of two free parameters (flux and polarized flux) . 
A Markov chain Monte-Carlo method was used to esti- 
mate constraints on the temperature profile as a function 
of radius. Specifically random steps were taken in the 
8-D parameter space and accepted if the proposed point 
was higher probability or rejected according to a proba- 
bility set by how less likely the proposed point was than 
the current point. After many samples, the list of points 
in the chain provides a sampling of the likelihood func- 
tion. For each point in the chain the temperature pro- 
file T(r) was constructed. The resulting distribution of 
temperature profiles is shown in Figure 1, where all pa- 
rameters have effectively been marginalized over. There 
is a dramatic reduction in the uncertainties if the cen- 
tral point source flux is assumed to be known both in 
terms of polarized emission and the SZ signal; this sim- 
ply reflects the fact that the polarized intensity is highly 
concentrated, and therefore can be partially mimicked by 
a point source. Given the importance of the central flux 
it seems prudent to simply marginalize over the unknown 
flux and accept that the temperature uncertainties will 
be on the order of <~ 1 keV, comparable to X-ray spec- 
troscopic temperature uncertainties. 

Real clusters are unlikely to follow ideal spherically- 
symmetric /3-profiles, and so we repeated the calcula- 
tion for clusters extracted from numerical simulations 
produced by the Hydra Consortium 8 and made publicly 

8 http://hydra.susx.ac.uk 
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Fig. 2. — Maps of thermal SZ (solid) and polarized inten- 
sity (dotted) for the most massive (left) and second most massive 
(right) clusters in the simulation archive. Outermost contours are 
1 fjiK at 30 GHz and they increase in factors of 10 inward. The 
polarized intensity assumes a central 1 Jy source, and that the 
conversion to temperature units is 1 /iK for 0.43 fijy per square 
arcminute. The direction of polarization at each point is tangential 
to the cluster center. 



available by the Virgo Consortium 9 on the web. The sim- 
ulation data shows deviations from the ideal /3-profilcs 
that are expected in nature . These simulations have been 
used (|da Silva et al. ll2000|) to study the statistics of SZ 
maps. Maps of two clusters, both in SZ and polarized 
emission, are shown in Figure 2. We repeated the above 
calculation, generating synthetic data, fitting to two j3- 
profiles and comparing to the true temperature profile. 
Results were qualitatively similar to those for the case 
of the spherical j3 clusters. In the simulated clusters the 
mass- weighted temperature within three arcminutes was 
recovered with statistical errors of roughly 1 keV. 

We have not attempted a detailed search to optimize 
the spatial sensitivity or frequency coverage. Contami- 
nating and confusing sources will be less or more impor- 
tant at different spatial scales or observing frequencies. 
At low frequencies there could be confusion with diffuse 
synchrotron emission by the cluster or extended jet struc- 
tures, while higher frequencies will run into both a rela- 
tive dearth of very bright central point sources and dust 
emission from our Galaxy as well a s distant star-forming 
galaxies ijFischer and Lange 111993^ . 

4. DISCUSSION 

We have shown that radio measurements alone can be 
used to reliably measure temperature profiles of the gas 
in X-ray clusters. The current generation of radio instru- 
ments is already sufficient in terms of sensitivity, but the 
presence of a 1 Jy source in the center of the field will 
no doubt make sub-mJy polarized studies challenging. 
Wit h the VLA per forming; at dynamic ranges exceeding 
10 5 llPerlev 1 119 99) . it is clear that the presence of the 
central source does not pose an inherent limitation to 
our method. 

With typical scattered fluxes of ~ 1% of the central 
source flux, it can be expected that the polarized scat- 
tered emission will exceed the /iK level for all central 
sources brighter than about 1 mJy at 30 GHz. This is 
brighter than the peak emission due to the cluster scat- 
tering of the primary CMB quadrupole, which is itself 
significantly larger tha n any anisotropics due to the rmal 
or kinetic SZ effects ijSazonov and Sunvaev Ifl999f) . A 
central point source is by far the single largest source of 

9 http:/ /virgo. Sussex. ac.uk/clusdata. html 



a quadrupole anisotropy in the radiation field for many 
clusters at radio frequencies. Therefore, detection of the 
scattered radio halo could provide an intermediate target 
for an attempt to measure these smaller signals. To the 
extent that the radio spectrum is constant over the age of 
the source, the spectrum of the polarized emission should 
be the same as that of the central source, providing an 
important consistency check. 

The nearly circular polarization pattern should pro- 
vide a simple template for improving the detectability 
of the signal, and the compact angular scale implies 
that the CMB anisotropy is not a serious contaminant. 
Many clusters show signs of r adio halos or radio relics 
ijGiovannini and Feretti 1 12002|) . These generally have 
steep spectra, so it may be difficult to observe the polar- 
ized emission reflected from the central source at wave- 
lengths longer than a few cm. Jets associated with the 
central source may be significant in the inner arcseconds. 
The largest contaminant will likely be the systematic 
problem of removing the very bright central source suffi- 
ciently accurately and stably so that it would not affect 
the measurements in the entire image. 

Nearby clusters are too large (in angular size) to have 
a significant surface brightness so the studies considered 
here are probably most suitable for clusters at redshifts 
z > 0.2. For studies of the SZ effect it should be noted 
that all sources above ~ 1 mJy will contaminate the im- 
mediate vicinity of the point source at /zK levels and 
therefore point source removal must be handled care- 
fully. If bright radio sources are significantly beamed, 
then some radio SZ observations could be contaminated 
by this reflected emission even for sources that do not 
appear to be very luminous. 

While we have outlined the potential for measur- 
ing temperature profiles, an equally exciting application 
would be to measure the time evolution of the radio flux 
from the central source. With cluster radii of roughly 1 
Mpc the light travel times are a few million years. With 
fiK sensitivities to polarized emission it will be possible 
to measure the age of the central source by simply finding 
an edge to the polarized emission; with independent tem- 
perature measurements, either from SZ spectral studies 
or X-ray spectroscopy, it would be possible to measure 
the time evolution of the central radio source. Surfaces of 
scattering corresponding t o constant central source emis - 
sion time are paraboloids ijSazonov and Sunvaev 1 12003) . 
and so a good mapping of the three-dimensional electron 
distribution will be necessary. 

In summary, the extremely high sensitivity and stabil- 
ity of current and near-future instruments should enable 
the use of the intracluster medium as a reflection neb- 
ula around bright point sources. This should provide 
some insight about the geometry of merging subclumps 
and a measure of the cluster electron density distribu- 
tion in addition to the important applications outlined 
above. The observations present some unique challenges, 
but the potential scientific rewards are compelling. 
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